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1 Introduction 

The formulation of relativistic two-body bound state wave equations and their 
relationship to quantum field theory began with the work by Eddington and 
Gaunt in 1928 1 . However, the large variety of approaches attempted in re- 
cent years shows that this problem still has no generally agreed-upon solution. 
Perhaps for this reason, most recent field theory books have skirted this topic. 
In his recent text, Steven Weinberg states "It must be said that the the- 
ory of relativistic effects and radiative corrections in bound states is not yet in 
satisfactory shape." 

Of course, this topic is often presented as covered by the manifestly covari- 
ant Bethe-Salpeter equation obtained directly from relativistic quantum field 
theory. Over the years, however, many problems have turned up to impede its 
direct implementation, mostly related to the central role played in it by the 
relative time or energy 0. These difficulties have led many authors to attempt 
reformulations. 

We describe here a recent approach resulting from " Two-Body Dirac Equa- 
tions" (emerging from Dirac's Relativistic Constraint Dynamics) that does sat- 
isfy many of the requirements one would demand of a treatment of the rela- 
tivistic two-body problem. We use its applications to QED bound states such 
as positronium, QCD quarkonia, and the nucleon-nucleon scattering problem to 
demonstrate the advantages of the approach. 



2 The One-Body Dirac Equation 

For a single body, the original Dirac equation for a single spin-one-half particle 
(nearly universally accepted) provides a successful bound state equation. The 
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free Dirac equation (7 • p + m)tp = serves as a relativistic version of New- 
ton's 1st law. When the four-vector substitution for electromagnetic interaction 
p^ — > p^ — Ap, and the minimal mass substitution for scalar interaction m — > 
m+5 are performed on it, one obtains a relativistic version of Newton's 2nd law 

(7 • (p - A) + m + S)*/j = 0. (1) 

The two-body Dirac equations of constraint dynamics successfully extend this 
one-body minimal coupling form to the interacting two-body system. 

3 Two-Body Dirac Equations from Constraint 
Dynamics 

In the 1970's, Todorov, Kalb and Van Alstine, and Komar independently used 
Dirac's constraint mechanics to attack the relativistic two-body problem for 
spinlcss particles 4 . By covariantly controlling the relative time variable these 
authors eliminated negative norm states as well as circumvented the no-interaction 
theorem of Currie, Jordan, and SudarshanJS] that had discouraged further work 
in this area since the 1960's . By combining constraint dynamics with par- 
ticle supersymmetries, Crater and Van Alstine extended those works to pairs 
of spin one half particles to obtain two-body quantum bound state equations 
that correct not only defects in the Breit equation but those in the ladder ap- 
proximation to the Bethe-Salpeter equation as well. [H|These Two-Body Dirac 
Equations of constraint dynamics possess a number of important features (some 
of which arc unique) which provide an alternative formulation of fundamental 
ficld-thcorctic results (while yielding standard perturbative spectra) and correct 
defects in phcnomcnological applications that result from patchwork introduc- 
tion of interactions : They a) provide a three dimensional but covariant re- 
arrangement of the Bethe-Salpeter equation, b) yield simple three-dimensional 
Schrodinger-like forms similar to their nonrelativistic counterparts, c) contain 
spin dependences determined naturally by their incorporation of Dirac's one- 
body structures, d) contain well defined strong potential structures that pass 
the necessary test that they reproduce correct QED perturbative results when 
solved nonperturbatively, e) in phcnomcnological applications make unneces- 
sary the ad hoc introduction of cutoff parameters generally used to avoid singu- 
lar potentials and f) have relativistic potentials which may be related directly 
to the interactions of perturbative quantum field theory or (e.g. for QCD) 
may be introduced scmiphcnomcnologically 7 .These equations provide a non- 
perturbative or strong-potential framework for extrapolating perturbative field 
theoretic results into the highly relativistic regime of bound light particles in a 
quantum mechanically well defined way. 

3.1 World Vector and Scalar Interactions 

The constraint formalism is embodied in a system of two coupled, compatible 
Dirac equations on a single wave-function. For particles interacting through 
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world vector and scalar interactions the two-body Dirac equations take the 
general minimal-coupling form 

Slip ee 751 (7 1 -(p 1 -ii)+mi + 5i)V = 

S 2 i> = 752(72 ' (P2 - M) + m 2 + S 2 )t[) = 0. (2) 

The two equations are compatible in the sense that, 

[S 1) S 3 ]il> = 0. (3) 

This condition is satisfied as a result of the presence in these equations of 
spin supersymmetries, a relativistic 3rd law, and covariant restrictions on the 
relative time and energy. Its direct dynamical consequence is the automatic 
incorporation of correct spin-dependent recoil terms . 

For the case where only vector interactions are present, the compatibility 
condition is most naturally satisfied by vector potentials in the "hyperbolic" 
momentum and spin-dependent forms 

ii = [1 - cosh(0)]pi + sinh(g)j5 2 -UdexpQ- 72 ) 72 (4) 

A 2 = [1 - cosh(£)]p 2 + sinh(£/)pi + - (dexpG ■ (5) 

In this case, the compatibility condition enforces a "third law" with the two 
constituent potentials actually depending on only one invariant .4. through the 
interaction function Q {A) . The explicit form of Q (A) follows both from com- 
parison with classical Wheeler-Feynman electrodynamics or with QED through 
leading term summation of ladder, cross-ladder and constraint diagrams 8 

1 2A 

Q{A) — — -log(l ) ; w the total cm. energy . (6) 

Compatibility also requires that the relative time be covariantly controlled 
through interactions depending only on x±, a covariant spacelike particle sepa- 
ration variable perpendicular to the total momentum P 

A =A(x ± ) 

p 

xf = x^ + P^fP ■ x), P ee — is a time-like unit vector. (7) 

w 

For lowest order electrodynamics, 

A = A(x ± ) = :r=Jx 2 ,. (8) 

r v ^ 

For quark-models, we must include scalar potentials Si . When appearing 
with vector interactions they depend not only on two invariant mass potential 
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functions M\(x±), M 2 (x±), related to each other through one invariant function 
L(x±) but also on the vector interaction through Q(A(x±)) 



Si 



Mi - mi 



-expQ(A)j 2 • 



S 2 



M 2 -m 2 + -cxpG(A)j 1 ■ 



M 2 ' 
dM 2 



(9) 



Mi ' 



Ml - Ml 




Mi = mi cosh L +m 2 sinh L 
M 2 = m 2 cosh L + mi sinh L 



;3RD LAW. (10) 



The counterpart to the invariant A for scalar interactions is S with the form of 
L = L(S(xj_),A(x±)) with m w — m\m 2 jw from 



Rctardative effects are already included through the cm. energy dependences of 
the potential structures. Although the potential forms in these equations may 
seem unfamiliar, expansion of the resulting classical dynamics in 1/c around 
the nonrelativistic limit shows that it is canonically equivalent to order 1/c 2 to 
the dynamics generated by the corresponding single-quantum exchange in field 
theory. 

3.2 Manifest Covariance and Well-Defined Quantum-Mechanical 
Behavior 

For our applications we use a Pauli reduction to bring our equations to the 
covariant Schrodinger-like form (with p the relative momentum) 



incorporating exact two-body relativistic kinematics through the eigenvalue 
b 2 (w) = (w 4 - 2(m 2 + ml)w 2 + (mf - m\^) 2 ) j Aw 2 = e 2 - m 2 w in terms of m w 
and e w = (w 2 — m 2 — m^)j2w respectively the mass and energy of the fictitious 
particle of relative motion. This Schrodinger-like equation is not only mani- 
festly covariant but quantum mechanically well defined: one can solve it non- 
perturbatively in both QED and QCD bound state cases for which every term 
in the quasipotential $ w (ai, <r 2 ,p±, A(r), S(r)) is less singular than — l/4r 2 (in 
contrast to all reductions of the Breit equation and many reductions of the 
Bethe-Salpeter equation). It involves at most 2 coupled wave equations but 
all portions of the 16 component wave function play essential roles in spectral 
calculations, cither directly or through the strong potential structures that they 
generate when they are eliminated. The explicit forms of the spin dependent 
potentials that appear in the Pauli-form quasipotential Q w are dictated by the 
interaction structure of the original two coupled Dirac equations and are not 
put in by hand. 



M 2 = m 2 +exp^(^)(2m tu S' + 5' 2 ) ;i 



= 1,2. 



(11) 



(p 2 +$ w (ai,<T2,P±,A(r),S(r)))il> = b 2 {w)ip 



(12) 
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The Schrodinger-like form of the Two-Body Dirac equations takes the min- 
imal coupling form 

(p 2 + (m, + S) 2 - (e w - A) 2 + 3>s P (vi,<r2,P±,A(r), S(r)))V = (13) 

in which 

$ S p = ®Dir -P+ $D2 + <S>SoL ■ ((7l + cr 2 ) + <S>SOdL • {a X - (T 2 ) + $SPO^ ' (0-1 X CT 2 ) 
-t-^SSfT! • (72 + $Tffi ' r<72 ■ f + $ DT? ' P "! ' ^02 • f. (14) 

The $i = <i>i (A, 5, w) are not independent but are all determined in terms of 
A, S through the Pauli reduction. 

We have checked analytically and numerically that our strong potential terms 
do not lead to spurious results by solving them nonperturbatively to obtain 
agreement with the standard fine and hyperfine spectra of perturbative QED. 
For example, for the singlet positronium system with A = —a/r our fully 
coupled system of 16-component equations S±ip = S21I) = is exactly solvable 
with total energy|5] 

w = m (2+2/(l+« 2 /(n+((;+l/2) 2 - a 2 ) 1 / 2 -;-l/2) 2 ) 1 / 2 ) 1 / 2 =m(2-^-^) ground statc 

(15) 

Such validation should be required of all candidate equations for nonpertur- 
bative quark model calculations and other semiphenomenological applications 
when their quark-model kernals are replaced by ones appropriate for QED. 
Otherwise, how can one trust the short distance spectral contributions obtained 
when applied in QCD? No other approaches have yet fully passed this test. In 
fact Sommerer et al 1101 have shown that all of the well-known quasipotential ap- 
proaches (e.g. the Blankenbecler- Sugar equation, the Gross, and Kadeshevsky 
approaches) fail this crucial test even for the ground state. However, by varying 
parameters permitted by the non- uniqueness of the quasipotential approaches, 
Sommerer et al obtained a quasipotential model that does reproduce this ground 
state numerically. We have successfully extended the check in our equations and 
found numerical agreement with the perturbative results for a range of angular 
momentum and radial states and for unequal masses. [H]. 

The invariant forms for A or S in our equations follow both from relativis- 
tic classical field theory and from the perturbative treatment of corresponding 
quantum field theories[y] through Sazdjian's derivation of the constraint equa- 
tions as a "quantum mechanical transform of the Bcthe-Salpeter equation" |11|. 

4 Two-Body Dirac Equations in Meson Spec- 
troscopy 

4.1 The Adler-Piran Potential 

We obtain a constraint version of the naive quark model for mesons from a 
covariant adaptation of a static quark potential due to Adler and Piran. These 
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authors use an effective non-linear field theory derived from QCD to find 12 

V AP (r) = A(£/(Ar) + U ) (=A + S). (16) 

The original Vap is nonrelativistic, appearing as the sum of world vector and 
scalar potentials in the nonrelativistic limit. Vap incorporates a running cou- 
pling constant form in coordinate space at short distance A[/(Ar << 1) ~ 
1/rlnAr and includes linear confinement plus subdominant logarithm terms at 
long distance 

V A p(r) = A(ciAr + c 2 log(Ar) + -p= + ^- +c 5 ), Ar > 2. (17) 

VAr Ar 

When the nonrelativistic quark model is constructed with realistic poten- 
tials such as the Adler-Piran potential or Richardson potential it fails for light 
mesons. Beyond a certain limit, it gives meson masses that increase with de- 
creasing quark mass. However, this effect is entirely absent from the light 
meson spectra produced by the Two-Body Dirac Equations with their relativis- 
tic kinematics and QCD-determined relativistic potentials.. 

4.2 Relativistic Naive Quark Model 

We reinterpret the static Vap covariantly by a) replacing the nonrelativistic r by 
\fx^_ and b) parcelling out the static potential Vap into the invariant functions 
A{r) and S(r) as follows 

A = exp(~(3r)[V A P --} + - + —, S = V AP + — - A. (18) 
r r r r 

This partially phenomenological step ensures that at short distance the potential 
is strictly vector while at long distance the vector portion is strictly Coulombic 
with the confining portion at long distance (including subdominant portions) 
strictly scalar. Note that in our equations, once A and S have been determined, 
so are all the accompanying spin-dependent interactions. Here our approach 
is that of a naive quark model since we ignore flavor mixing and the effects of 
decays on the bound state energies. 

4.3 Meson Spectroscopy 

Our results are spectrally quite accurate, from the heaviest upsilonium states 
to the pion. Notable exceptions are light meson orbital and radial excita- 
tions and their spin-orbit splittings. We include only selected portions of the 
whole table here (see hep-ph-0208186 for the complete results). The param- 
eter values we obtain from the best fit are nib = 4.877, m c = 1.507, m s = 
0.253, m u = 0.0547, m d = 0.0580 GeV along with A = 0.216, AU = 1.865 
GeV and f3 — 1.936. We find that the static Adler-Piran potential, having 
a close connection with quantum chromodynamics (QCD), gives a good fit to 
the mostly nonrelativistic bottomonium spectrum. We obtain 9.453, 9.842, 
9.889,9.921,10.022 GeV for the ground state and first orbital and radial states 
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of 9.460,9.860,9.892,9.913,10.023 GeV. For the B mesons the fit results from 
the fact that for them our equations essentially reduce to the one-body Dirac 
equation We obtain for the 5(5.279), B*(5.325), B s (5.369), B* (5.416) mesons 
the results 5.273,5.321,5.368,5.427 GeV. The goodness of the fit to the charmed 
quark mesons shows that the equations perform well in the semirelativistic re- 
gion. The ground state ( 1 S , , 3 S'i) and first orbital ( x Pi, 3 Pq, 3 Pi, 3 P2) and ra- 
dial (2 1 S ,2 3 S 1 ) excitations of 2.980, 3.097, 3.526, 3.415,3.510,3.556,3.594,3.686 
GeV match our theoretical values of 2.978, 3.129,3.520, 5.407,3.507, 3.549,3.610,3.688 
GeV quite nicely. The same thing happens for the £>(1.865), D*(2.007), D s (1.968), £>*(2.112) 
mesons for which we obtain the values of 1.866,2.000,1.976,2.123 GeV. The good- 
ness of the accompanying fit to the lighter mesons (with the same two invariant 
potential functions used for the entire spectrum) is due to exact two-body rela- 
tivistic kinematics combined with the minimal interaction and strong potential 
structures of our equations for vector and scalar potentials. For example, 
for the ground states AT (494), if*(892), </>(1.019), 7r(140),and p(767), we obtain 
0.492,0.910,1.033, 0.144, 0.792 GeV. The results for the light meson orbital and 
radial excitations and their spin-orbit splittings are mixed. For example the 
6i(l. 231), o (1.450),ai(l. 230), a 2 (l. 318) meson fits of 1.392,1.491,1.568,1.310 are 
quite uneven and the radially excited 7r(1.300) and p(1.465) results of 1.536 
and 1.775 GeV are quite far off the mark. On the other hand the results 
1.319,1.533,1.493 GeV for the cf) orbital excitations are reasonable. In the fu- 
ture, we plan to use a coupled channel formalism to investigate the origin of 
some of these problems. 

The strength of the Two-Body Dirac approach is that using it , with just 
two parametric functions A and S we are able to obtain an overall fit about 
as good as that obtained by Godfrey and Isgur|13|. who used six parametric 
functions, basically one for each type of spin dependence. 

As a bonus, we find that the pion is a Goldstone boson in the sense that 
m v (m q — > 0) — > while the p and excited 7r have finite mass in this limit (see 
hep-ph/0208186 for complete plots). We have shown elsewhere that (0|c? m j^ |7r) = 
(Mi(0) + M 2 (0)) = (mi + rri2)Tr('y 5 ip) supporting our claim that our potential 
model incorporates aspects of a spontaneous breakdown of chiral symmetry. 

5 Two-Body Dirac Equations in Nucleon-Nucleon 
Scattering 

The two-body Dirac equations of constraint dynamics provide a natural method 
for extending nonrelativistic phenomenological treatments to the relativistic 
domain with effective potentials determined from the standard field-theoretic 
treatment of meson exchanges. This provides a severe test for the strong 
potential terms that turned out to be essential for the treatment of QED 
and QCD bound states. The mesons we include are the pseudoscalar mesons 
?r(135), 7?(548), 7/(952) the vector mesons p(770), w(776), 0(1020) and the scalar 
mesons cr(600), a (980), /o(983). The it, p, and the ao are isovector mesons while 
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the rest are isoscalar mesons. To incorporate these nine meson exchange forces 
into the two-body Dirac equations we need to generalize the interactions con- 
tained in our equations to include others beyond world-scalar and world- vector. 

5.1 Two-Body Dirac Equations for General Covariant In- 
teractions: The role of supersymmetry 

The detailed forms of interaction in our equations are actually the consequences 
of supersymmetries in our interacting two-body system. As an illustration, we 
review the derivation of Eq. J2J) . Define theta matrices in terms of Dirac matrices 
M = i072~7 5 7f [x = 0, 1,2,3, 6I5 = In the "correspondence" limit 

in which the 0's become Grassmann variables, the Dirac equation becomes a 
constraint imposed on both bosonic (j>) and fermionic (9, 9$) variables: 

S ip = (p-9 + m6 5 )ip = 0;=>S„ = (p-9 + m9 5 ) « 0. (19) 

For a single free particle, this "classical Dirac-Equation" constraint is supersymmetric|H], 
under the supersymmetry generated by p ■ 9 + \] —p 1 9^ which however does not 
leave the position four-vector x invariant . However, it does leave the "zit- 
terbewegungless" position variable x^ = x^ + i9^9^jm invariant. In the 
presence of scalar interaction M = m + S , i M becomes the "self-referent" 
form x^ — x^ + i9^9$jM{x) . The "pseudoclassical" Dirac dynamics is then 
governed by the supersymmetric system of constraints 

S = p ■ 9 + M(x)9 5 « 0, -{S, S} = H=p 2 + M 2 {Z) w 0. (20) 

i 

Since 9\ = , the expansion of the self-referent form truncates so that M(x) = 
M(x)+idM(x)-99s/M(x). Upon quantization, Eas. (|20|) then turn into the Dirac 
Equation and its standard square when the Grassmann variables become theta 
matrices while dynamical variables x and p become their operator counterparts. 
Thus the supersymmetry that preserves J is a natural feature of both the single- 
particle Dirac equation for the free case and its standard form for external scalar 
interaction. 

For two particles we introduce interactions that preserve such a supersym- 
metry for each spinning particle through the replacement 

m« -> Mi(xi - x 2 ) -> Mi(xi - i 2 ) = Mi, i = 1,2. (21) 

The Grassmann Taylor expansions of the Mi truncate leading to 

Sxi) = (9 1 -p + e 1 9 1 -P + M 1 9 51 -idL-9 2 9 52 9 51 W = 0, 

S 2 *P - (-02 ■p + e2&2-P + M 2 52 +idL- 01052051)^ = 0, (22) 

with the invariants Mj and L(xj_) related by Eci. ((TU)l . Eq.® becomes Ea. (|2*2l 
when restricted to scalar interactions. The consequences of pseudoclassical su- 
persymmetries are the extra spin-dependent recoil corrections to the ordinary 
one-body Dirac equations, essential for the compatibility of the two equations. 
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5.2 Hyperbolic Form of the Two-Body Dirac Equations 
for General Covariant Interactions 

We introduce general interactions by recasting the minimal interaction forms 
of the two-body Dirac equations into ones that generalize the hyperbolic forms 
we encountered in the treatment of scalar interaction. In the scalar case, if we 
begin with the two constraints 

SiV> = (Sio cosh(A)+S 20 sinh(A))V> = 0, S 2 ^ = {S 20 cosh(A)+5i sinh(A))^ = 
we find that Eas. (|22|) are equivalent to the linear combinations 

Sitp = [cosh(A)Si + sinh(A)S 2 ]^ = 0, S 2 ^ = [cosh(A)S 2 + sinh(A)Si]V> = 

(24) 

in which the interaction appears through the invariant matrix function A = 
—95i9§ 2 L(x±_) while the Sto — (pi ■ 9i + ™#5i) are free Dirac operators. The Si 
constraints (and hence the Si ) are a compatible pair for general A: [Si, S 2 ]?/> = 
(and [iSi,5 2 ] , = 0) provided only that A = A(x±) . Consider the four polar 
and four axial interactions. For the polar interactions we find the forms A(a;j_) = 
— L(x±)85i952 for scalar, J{x±)P ■ Q\P ■ 9 2 and G(x±)9i± ■ 9 2 ± for time-like and 
space-like vector respectively and A(xj_) = J-(xj_)6i± ■ 2 ±05i6*5 2 P • 9\P ■ 9 2 for 
polar tensor. The constraint equations for vector and scalar interactions given 
in Eq.Q come from L and the Feynman gauge combination Q = — J when 
T =0. The Si constraints for the axial counterparts are defined as above but 
the Si linear rearrangements have a minus sign in place of a plus sign. For 
them the invariant A forms are A(xj_) = C(x±)/2 for pseudoscalar , A(xj_) = 
H{x±)P ■ 9\P ■ 9 2 9^\9^ 2 and I(x±)9i± ■ 9 2 ±95i05 2 A(x±) for time-like and space- 
like pseudovector respectively and A(x_l) = Y(x±)9i± ■ 9 2 ±P ■ 9\P ■ 9 2 for axial 
tensor. Some of these may find application in meson spectroscopy in addition 
to the usual scalar and vector interactions. 

5.3 Strong Potential Forms for Nucleon-Nucleon Scatter- 
ing 

To represent meson exchanges, we require pseudoscalar interactions as well as 
vector and scalar. Thus we must use 

A(x±_) = -L(x x )9 5 i9 52 + 0(x_l)0i • 9 2 - C(x x )/2. (25) 

with the electromagnetic four- vector (Feynman gauge) condition J(x±) = —G(x±) 
relating time and space-like components {9\-9 2 = 9i±-9 2 ±—P-9iP-9 2 ).Reduction 
of the coupled Dirac equations (|24() to Schrodinger-like form for these combined 
interactions produces |14j for equal masses 

$w — » $si + $D + (§so + $SOTVi-ra 2 -f)L-(a 1 +a 2 ) + $ssVi-o'2 + ®TO-i-fa 2 -f. 

(26) 
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in which L, C, Q fix all separate quasipotential pieces. Note that the quadratic 
nature of many of the strong potential terms (e.g. S 2 , A 2 , L' 2 , Q' 2 , L'C etc.) 
could lead to disastrous results for large coupling-constants. 

Beyond the limitation to the above nine mesons, one must choose how the 
corresponding nine Yukawa potentials are included in the three invariant func- 
tions C,L,Q = — J . We structure our strong potential terms by assuming that 
in Q = -\ log(l - ^A) we take A = A for A < and 

. w ,An. 
A = — arctan( ) , A > 0. * * 

7T W 

2 - - exp(-TOpf) 2 expj-TOj^F) 2 exp(-m f) 
A = g p T 1 -T 2 + g w + 9(j> (27) 

For the invariant L we use Ea. (|llf) for S > and 
1 2 9 

5 = - g r xp( ; mgf:) - ^ exp( ~ m/ ° f) - fli^i ■ r 2 exp( 7- r " ) . (28) 

The modifications (**) of our strong potential terms for large repulsive vector 
and large attractive scalar interactions lead to corresponding changes in the 
quasipotential portions. For the pseudoscalar invariant function C we use 

1 2 ^ „ cxp(-TO 7r f) 2 cxp(-m t) f) 2 exp(-m ?? /f) 

k = — SVT1-T2 hg^ Vg^, = J. (29) 

w r r r 



We model effects of form factors by replacing r {=\J%\) by f = \Jr 2 + r 2 . In 
addition we take into account that the vector mesons may have an anomalous 
"magnetic moment" type of coupling. The net effect is to include pairs of 
additional vector and scalar Yukawa interactions but with opposite signs. 

For 3 Si n — p and 1 Sq n — p scattering we obtain excellent phase shift fits 
in the energy range from 1 to 350 MeV . Examination of some of the other 
scattering states shows, however, that the model needs improvement through in- 
clusion of a) world tensor coupling, b) pseudovector coupling of the pseudoscalar 
mesons and c) the off mass shell effects of the vector meson couplings. 



6 Concluding Remarks 

Nonperturbative solution of our two-body Dirac equations (with all of their 
strong-potential structure) in QED has demonstrated that these relativistic 
wave equations reproduce the field-theoretic perturbative spectral results thereby 
increasing our confidence in their use whenever Coulomb-like potentials play 
a significant role in dynamics in field-theoretic or phenomenological applica- 
tion. They have been successfully applied in QCD with spectral results as good 
as those of the most popular approach but using just two invariant potential 
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functions (with Goldstone boson behavior as a bonus). Their application to 
NN— scattering necessitates their extension to include general covariant interac- 
tions (achieved through a natural hyperbolic structure present in them). Phase 
shift results obtained from them look promising. In the future, the covariant 
and local Schrodinger-like structure of the equations that make them simple 
to implement may allow us to combine them with formalisms originally devel- 
oped for the nonrelativistic Schrodinger equation e.g. the microscopic theory of 
meson-meson scattering |15] and the unitarized quark model 16 . 
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